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Abstract. Assuming the Riemann Hypothesis, we show that infinitely often 
consecutive non-trivial zeros of the Riemann zeta-function differ by at most 
0.5155 times the average spacing and infinitely often they differ by at least 
2.69 times the average spacing. 



1. Introduction 

Let C,{s) be the Riemann zeta-function. Assuming the Riemann Hypothesis, we 
denote the non-trivial zeros of ({s) as p = ^ ± ij where 7 e K. It is weh known 
that, for T > 10, 

0<7<T 

Hence, if we let 7 < 7' denote consecutive ordinates of the zeros of C(s), we see 
that the average size of 7' — 7 is 27r/ log 7. Normalizing, we let 

(7' - 7) log 7 



A := limsup 

and 



7>o 27r 



. „ (7'-7)log7 
u := iimmi 

and we observe that /i < 1 < A. It is expected that there are arbitrarily large 
and arbitrarily small (normalized) gaps between consecutive zeros of the Riemann 
zeta-function on the critical line; in other words, that fi ~ and A = +00. In this 
note, we prove the following theorem. 

Theorem 1.1. Assume the Riemann Hypothesis. Then A > 2.69 and fi < 0.5155. 

We briefly describe the history of the problem. Very little is known uncondi- 
tionally; however, Selberg (unpublished, but announced in [11]) has shown that 
fj. < 1 < X. Assuming the Riemann Hypothesis, numerous authors [3 El [6l [7l [9] 
have obtained explicit bounds for and A. Theorem 11.11 improves the previously 
best known results under this assumption which were /i < 0.5172 due to Conrey, 
Ghosh & Gonek [2] and A > 2.63 due to R. R. Hall |;5J. Assuming the generalized 
Riemann Hypothesis for the zeros of Dirichlet L-functions, Conrey, Ghosh & Gonek 
[3] have shown that A > 2.68. Their method can be modified (see [101 and [1]) to 
show that A > 3. 



2000 Mathematics Subject Classification. Primary 11M26; Secondary 11M06. 

HMB is supported by an EPSRC Postdoctoral Fellowship. MBM is supported in part by a 
University of Mississippi College of Liberal Arts summer research grant. NN is supported in part 
by a an NSERC Discovery grant. 



2 



H. M. BUI, M. B. MILINOVICH, AND N. NG 



Understanding the distribution of the zeros of the zeta-function is important for a 
number of reasons. One reason, in particular, is the connection between the spacing 
of the zeros of C(s) and the class number problem for imaginary quadratic fields. 
This is described by Conrey & Iwaniec in [4]; see also Montgomery & Weinberger 
[8j . Studying this connection led Montgomery [6] to investigate the pair correlation 
of the ordinates of the zeros of the zeta-function. He conjectured that, for any fixed 
< a < /3, 



^ ' sin TTit ^ 2 



0<7,7'<T 



du. 



Here 7 and 7' run over two distinct sets of ordinates of the non-trivial zeros of 
C(s). Clearly, Montgomery's conjecture implies that = 0. Moreover, F. J. Dyson 
observed that the eigenvalues of large, random complex Hermitian or unitary matri- 
ces have the same pair correlation function. This observation (among other things) 
has led to a stronger conjecture that the zeros of the zeta-function should behave, 
asymptotically, like the eigenvalues of large random matrices from the Gaussian 
Unitary Ensemble. These ideas lead to the conjecture that A = +00. 

2. Montgomery & Odlyzko's method for exhibiting irregularity in 
the gaps between consecutive zeros of 

Throughout the remainder of this note, we assume the truth of the Riemann 
Hypothesis. 

Let T be large and put K = r(logT)-2. Further, let 



where 



3c 



2sin(7rci^) 



Trlogn 

and A(-) is von Mangoldt's function defined by A(n) = logp if n ~ for a prime 
p and by A(n) = 0, otherwise. In |7], by an argument using the Guinand-Weil 
explicit formula for the zeros of C(s), Montgomery & Odlyzko show that if /i(c) < 1 
for some choice of c > and a sequence {an} then A > c and that if h{c) > 1 for a 
choice of c > and a sequence {an} then pL < c. In particular, for any such choices 
of c and {an}, their method proves the existence of a pair of consecutive zeros of 
C(s) with ordinates 7 < 7' in the interval \T/2, 2T] which satisfy 7' — 7 > and 

7' - 7 < liff^, respectively. 

Conrey, Ghosh & Gonek [2] have given an alternative, and much simpler, way 
of viewing this problemQ Let 



A{t) = E ^kk- 



k<K 



^Mueller [9] was the first to observe that this method gave a lower bound for A. In [2], it was 
noticed that the method also gave an upper bound for and was equivalent to the method in [7]. 
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be a Dirichlet polynomial and set 

Mi^ \A{t)\^dt and IVhic) ^ / V \A{-f + a)\^ da. 

Jt J-ttc/ log Tj,<^<2T 

Then, clearly, M2(c) is nionotonically increasing and M2(/i) < Mi < M2{X). There- 
fore, if it can be shown that M2(c) < Mi for some choice of A{t) and c, then A > c. 
Similarly, if M2(c) > Mi for some choice of A{t) and c, then jjL < c. Using standard 
techniques to estimate Mi and M2(c), it can be shown that 

M2{c)/Mi^h{c)+o{l). 

Hence, this argument is seen to be equivalent to Montgomery & Odlyzko's method, 
described above. Moreover, we note that this formulation of the method suggests 
that we should choose a test function A{t) which is small near the zeros of C(s) to 
exhibit large gaps between the zeros of the zeta-function and a test function A{t) 
which is large near the zeros of C(s) to exhibit small gaps. 
In [7j, Montgomery & Odlyzko make the choices of 



-L/(^) and = 



where / is a continuous function of bounded variation on [0, 1] normalized so that 
/o l/P = 1 ^ii'i ^{k), the Liouville function, equals (— l)^'*^'; here, ^{k) denotes 
the total number of primes dividing k. By optimizing over such functions /, the 
values /I < 0.5179 and A > 1.97 are obtained. The authors of [7] choose / to be a 
certain modified Bessel function and they mention this is close to an optimal choice. 
In [2], Conrey, Ghosh & Gonek choose the coefficients 

dr{k) , X{k)dJk) 

"^■^^ "'^^^^ 
where dr{k) is a multiplicative function defined on integral powers of a prime p by 

, _ r(fc + r) 

'■^^ ^ " T{T)k\ ■ 

In this context, this becomes an optimization in the variable r. The choice r = 1.1 
yields ^ < 0.5172 and the choice r = 2.2 yields A > 2.337. 

In order to prove Theorem II. li we combine the approaches of [7j and [2] . We 
choose the coefficients 

dr{k) , log K/k\ , \{k)dr{k) ,\ogK/k\ 

for sufficiently smooth functions /. This variant allows us to optimize over both r 
and /, rather than over just r or /. 

We now provide further insight into the choice of these coefficients. For simplic- 
ity, suppose / is a polynomial. Since, for Re s > 1, 

k>l k>l ' 

with our choice of coefficients we see that the test function A{t) approximates 
+ ity and C(l -|- 2itY / (^{j -\- ity , respectively, and should have the desired 
effect of making A{t) small (respectively large) near the zeros of ((s). Moreover, 
when we multiply dr{k) by / ( ) then A(t) behaves like a linear combination 
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of + ity and its derivatives and an analogous comment applies to the other 
case. The presence of the function / leads to improved numerical results for bounds 
for and A. 



3. Large Gaps: a lower bound for A 

In this section, we establish a lower bound for A by evaluating h{c) with the 
coefRcients 



f { log K/k -\ 



where / is a continuous, real- valued function of bounded variation on £^[0, 1]. In 
what follows, we assume that r > 1 so that dr{mn) < dr{m)dr{n) for m, n e N. It 
is well known that, for fixed r > 1, 

(3.1) J2 ^ - ^.(log:^)'^^ + OiilogTf-') 

k<x 

uniformly for x < T; here Ar is a certain arithmetical constant (the exact value 
is not important in our argument). By partial summation, recalling that K = 
T(logr)~^, we find that the denominator in the ratio of sums in the definition of 
h{c) is 

k<K k<x 

by p.ip . By the variable change u = 1 — 

J2 \akf^ Ay {log KY' [\l-uY'-'f{ufdu + Of.r{{\ogTf-') 

k<K 



= Ay{\ogTy / {i-uY-^f{ufdu + Of.r,e{{\ogTy-^+') 

Jo 

where e > is arbitrary. 

We now evaluate the numerator in the ratio of sums in the definition of h{c). If 
we let 

nk<K 

then a straightforward argument shows that 

2 dr{k)dr{kn)K{n) ^ no^) f (]2^]£hh) sin iirc'^^) 

\ ' ^ ZL^ Jq„ ^ J \ \ogK I ■> y \ogK ) \ log T / 

nk<K ° 

= ^ E ^-(^);^;(^^V (^)/(^)Bin(vrci^)+Q,,.((logT)---^) 

pk<K ^ 

^g^^ sin(.cj^) ^ '^fC-^)fC-2^)+0Mil0ETf-'). 

p<K ^ k<K/p 
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By Stieltjes integration and a variable change, the inner sum in the main term of 
the last expression for N{c) is 

r'/(^)/(^)<E^) 

^ k<x 

/K/p 1 

= Arr\\ogKff {l-uf-^f{u)f{u-^) du + Of,r{{\ogTf-'). 

By combining the above estimates and interchanging the order of summation and 
integration, we conclude that N{c) = M(c) + 0/^r((logr)'' where 

2Ar^ 2 2 ^ sin (ttc^) 

•^IS#K 2<p<K" ^ 

= — — (log^) / (l-"") /(") 2^ f\u-\^)du 

+ 0/,.,e((iogrr^-i+^). 

By the prime number theorem with remainder term, it follows that 

2^ /I" iogK>-J ^1 /I" log^^'^'^ + ^^^'^ViogrJ- 

By the variable change v = 1°^^ , the integral is 
sin (TTCwi^S^) /■« sinfTTCt;') 

7log2 V Jo V 

log ft: 

Hence, 

M(c) = (log (\l~uY"-^f{u) r^J^^^^f{u~v)dvdu 

71" Jo Jo V 

+ 0f,rAi^0gTf-'+^). 

Consequently, we find that 

^ ^ _ 2r j',{l^ur-'f{u)j;^^^f{u-v)dvdu 
(3.2) ^ ^ TT Sl{\-uY'-^f{uYdu 

+ 0/,.,e((logT)-i+^'). 

Choosing r = 3.1 and f{x) = l + 10a; + 39a;^, we obtain (by a numerical calculation) 
that /i(2.69) < 1 when T is sufficiently large. This provides the lower bound for A 
in Theorem ll.il 
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4. Small Gaps: an upper bound for ^ 

Since A(n)^ = 1 and X{pn) ~ —X{n) for every n G N and every prime p, we can 
evaluate h{c) using the coefficients 

_ X{k)dr{k) ^ f log K/k \ 

in a similar manner to the calculations in the previous section. Here, as above, 
/ e L^[0, 1] is a real- valued, continuous function of bounded variation. With this 
choice of coefficients, we obtain that 

(4.1) " TT J^il-ur-^f{uydu 

+ 0/,.,,((logT)-i+^'). 

Choosing r — 1.23 and f{x) = 1 + 0.99a; — 0A2x^, a numerical calculation implies 
that /i(0.5155) > 1. This provides the upper bound for /i in Theorem ll.il 



5. Some Concluding Remarks 

Theorem 11.11 offers the best known bounds for A and /i assuming the Riemann 
Hypothesis; however, we are still far from proving the conjectured values of /i = 
and A = oo. In fact, it known that this is not attainable using Montgomery and 
Odlyzko's method. Specifically, in [2], it is shown that h{c) < 1 if c < i and 
h{c) > 1 if c > 6.2. It would be interesting to better understand the limitations of 
this method and, in particular, if it can be used to show that /i < ^. 

We have not been able to prove that our bounds for A and jj, in Theorem 1 1.1 1 are 
the optimal bounds for our choice of coefficients 



dr{k)_fflogK/k\ , _ X{k)drik) ^f log K/k ^ 

-^n-i^iir) and flfc- ilTSiir; 



In the special case of r = 1, this optimization problem has been solved (in terms 
of prolate spheroidal wave functions). When r ^ I, the analogous optimization 
problem seems considerably more difficult. Instead of trying to solve it explicitly, we 
have instead chosen / to be a polynomial of low degree (< 6). Using Mathematica, 
we numerically optimize (|3.2p and (14. ip over the coefficients for each choice of c 
and r. 

Our numerical calculations seem to suggest that degree 2 polynomials work well; 
it does not seem like there is much to gain by increasing the degree of /. To 
demonstrate this phenomenon, we observe that one can recover the bounds for A 
and /i, in the case of r = 1, derived in |7] using degree 2 polynomials in place 
of prolate spheroidal wave functions (which are much more difficult to compute 
numerically). In this case, after some rearranging and a change of variables, our 
above calculations imply that 

h{cj = c± J ^ h o(l). 

Jo fi^Y du 

Letting f{x) = 1 + 0.46526a; — 0.46526a;^, a numerical calculation shows that 
/i(0.5179) > 1 and if we let f{x) = 1 + 17.9426a; - 17.9426a;2, then it can be 
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shown that ft. (1.97) < 1. These are essentially the (optimal) values obtained by 
Mongomery and Odlyzko in 7J when r ~ 1. 
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